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Abstract: Primary gamma rays emitted by extragalactic sources, such as blazars, will generate electromagnetic cascades 
in intergalactic space. These cascades proceed via electron-positron pair production and inverse Compton scattering 
on cosmic background radiation, mainly the cosmic microwave background (CMB) and extragalactic background light 
(EEL) fields. The existence of an extragalactic magnetic field (EGMF) could deflect electron-positron pair trajectories 
and scatter the cascade photons, possibly creating a halo around the source while suppressing the cascade flux collected 
by a detector We develop a semi-analytic model for the cascade process and apply it to combined GeV-TeV data on high- 
frequency-peaked BL Lacertae objects (HBLs) from the Fermi Large Area Telescope (LAT) and ground-based Cherenkov 
telescopes, comparing observation results with model predictions using a robust statistical framework. Lower limits with 
different confidence levels on the field strength of the EGMF derived from this procedure are discussed under various 
assumptions about the source livetime. 

Keywords: Gamma-Ray Astronomy, Extragalactic Magnetic Field, Electromagnetic Cascade, BL Lacertae Objects, 
Extragalactic Background Light, Fermi Gamma-Ray Space Telescope, VERITAS 



1 Introduction 

Extragalactic blazars emit gamma rays in both high-energy 
(HE, 100 MeV < E < 300 GeV) and very-high-energy 
(VHE, E > 100 GeV) bands. Due to the existence of 
the extragalactic background light (EBL), which spans over 
the optical to far-infrared wavelength range, gamma rays 
with energies above 10 GeV may be absorbed and pro- 
duce electron-positron pairs |r|. An electromagnetic cas- 
cade then develops via inverse Compton scattering of the 
pairs on the cosmic microwave background (CMB) and 
subsequent secondary pair productions. In the presence of 
an extragalactic magnetic field (EGMF), the charged pairs 
in the cascade will be deflected, spreading the cascade pho- 
tons in both the spatial and temporal distributions IJIO. 
The characteristic angular spread could create an apparent 
halo around the point source, and the time delay of cas- 
cade photons could appear in the observation of gamma- 
ray bursts or flaring blazars. Therefore, gamma-ray astron- 
omy of extragalactic sources provides a useful probe into 
the EGMF strength and configurations. 

According to previous studies f4^, '31 this method will be 
sensitive to EGMF below 10^^* Gauss, much lower than 
what any other measurement has achieved. For instance, 
from Faraday rotation measurements on extragalactic radio 
sources El?) or analysis of CMB anisotropy |8 9| only an 
upper limit ^ 10^^ Gauss on the field strength is obtained. 
The new EGMF window below 10^^^ Gauss is particularly 
interesting for the understanding of astrophysical magnetic 



fields. A primordial field within the window could be re- 
sponsible for generating the galactic and intra-cluster mag- 
netic fields |12|, and its own origin can be related to ei- 
ther the inflationary era or phase transitions in the early 
universe ifTJl . On the other hand, if the EGMF strength 
turns out to be zero, the astrophysical fields would have to 
be coming from seed fields produced locally via the Bier- 
mann battery mechanism or related processes llT4l . Hence 
a lower limit on the EGMF instead of a measurement could 
be already useful in clarifying the origin of all the magnetic 
fields we have in the universe today. 
One way for obtaining the lower limit, using gamma rays as 
a probe, would be to compare the cascade flux from VHE 
emission of a source with the actual HE observed spec- 
trum IITSi [161 . If the HE measured flux is lower than the 
zero-field cascade flux prediction, the EGMF would have 
to be non-zero to dilute the cascade photons into a spread- 
ing angle and thus suppress the collected flux. Deriving an 
EGMF lower limit in this way requires a realistic model 
of the electromagnetic cascade correlated with EGMF. Ex- 
isting simplified analytic models |17| and Monte Carlo 
simulations 1 1 8 1 set a lower bound for the field strength 
at 10^^® to 10^^^ Gauss assuming the studied sources to 
be active with unlimited livetime or 10~^^ to 10~^^ Gauss 
for the sources to be active for only ~ 3 years of simulta- 
neous HE and VHE observations 1(19112011211 . In this work 
we model the cascade semi-analytically and use our model 
predictions to place a lower limit on the EGMF using a sys- 
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Figure 1 : Sketch of the gamma-ray cascade geometry. 



tematic framework, with the blazar RGB J0710H-591 as an 
example for appHcation and data analysis. 

2 Model Description 

The geometry of the cascade is shown in Fig. [1] Primary 
gamma rays emitted by a blazar at distance L from the earth 
are absorbed after going through distance L' . The electron- 
positron pairs get deflected by the EGMF to angle Qd and 
upscatter CMB to secondary photons directed toward the 
detector at an incidence angle 0c- The emission angle of 
the primary photon at the source with respect to the line of 
sight is Bs ~ 9d — 0c- The difference in path length between 
the secondary photon and a direct photon that goes from the 
source to the observer in a straight line is 



AL = cAT = L' + ^Jl^ + L'^ - 2LL'cos9s - L (1) 

where AT is the time delay of the secondary photon and c 
is the speed of light. 

For a primary photon with energy e, we assume the elec- 
tron and positron produced in the pair-production process 
each carry e/2. Then the pair starts to lose energy contin- 
uously via inverse Compton scattering on the CMB while 
being deflected by the EGMF at the same time. An electron 
with Lorentz factor 7e on average scatters CMB photons to 
energy 47^60/8, where cq is the average CMB photon en- 
ergy ^ 6.4 X 10""* eV II22I . The energy loss process is then 
described by 



at 6 



(2) 



where ricMB is the CMB photon number density ^ 
411cm~'^ and ctt is the Thomson cross section ^ 6.65 x 
10~^^cm^. The Lorentz deflection process is, on the other 
hand 

d9d _ _c _ eB 

dt n JeTn-eC 

where ri = •je'meC^ / eB is the Larmor radius, me is the 
electron rest mass, and we have assumed the EGMF B 
to be perpendicular to the electron momentum. Eqs. |2] 
and [3] combine to give the deflection angle for an elec- 
tron/positron to go from Lorentz factor to 



eB 



< eo"-CMBO'T 



(7e ^ - 7eO^) 



(4) 



9d = arccos (sin^ 9 f cos 9dx) + cos^ 



(5) 



when the angle 9 / between B and the electron momentum 
is other than 7r/2. 

Combined with the geometry in Fig. [T]Eq. |5] uniquely de- 
termines 9s and 9c for a given set of B, 7eo, 7e, L', and 
9f, provided that 9c < Tr/2. The number of secondary 
photons between energies E and E + dE produced by the 
electron going from 7e + dje to 7e can be calculated from 
the CMB spectrum, replacing the CMB photon energy with 

3^^/(47?): 



dN{E,^c) 



cdtu' 



277tE^ 



dE 



874 /i3c372 (^e3E/i(^^kT) _ 

SlTTE'^nied'^c dE 



(6) 



where h is the Planck constant, k the Boltzmann constant, 
and T the CMB temperature at 2.73 K. Integrating over 
7e, 7eO (or equivalently e), L', and averaging over 9f gives 
the differential secondary photon flux for a certain field 
strength B as 



dN{E) _ 81ttE^7 



dE 16/i3ceoncMB 7 78 (e3S/(47|feT) _ 1) 



dje 



J d9fg{9f) J dej dL' 



A(6) 



-f{e,9s) 



exp (-^/L^ + L'^ - 2LL' cos9s/X{E) 



(7) 



where A(e) is the mean free path of a gamma-ray photon at 
energy e, depending on the specific EBL profile. g{9f) is 
the probability distribution of 6*/, which is sin 9 / for a ran- 
domly pointing field. f{e,6s) is the intrinsic spectrum of 
the source, and we integrate over e starting from 2^emeC^ 
to 200 TeV, as primary photons beyond that energy are 
mostly absorbed within 1 Mpc away from the blazar, and 
quickly deflected away by the relatively large magnetic 
field there. The upper limit on 7e is 100 TeV/(meC^) and 
a lower limit at lO'"^ is also placed as a practical matter for 
the numerical integration, since there is negligible CMB 
density beyond 3 meV and we are not interested in sec- 
ondary flux below 100 MeV. The integration limits on L' 
are enforced through observational cuts on AT and 9c via 
geometry shown in Fig.[T] 

The blazar intrinsic emission f{e,9s), the photon mean 
free path A(e), and the EGMF directional profile g{9f) are 
the inputs to the cascade model. In practice we model the 
blazar emission as boosted isotropic radiation | 



/o(l-/3cos0,) 



+ /o(l+/3cos6',)-"-'e-"e-'/^° (8) 



where the second term models a counter jet. To obtain a 
conservative prediction on the cascade flux we choose the 
optical depth profile of [241 which is relatively transparent 
for VHE gamma rays. g{9f) is taken as sin 6*/ as we have 
no prior assumption on the EGMF configuration. 
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Figure 2: Measured spectral points of RGB J0710+591 and 
predicted spectra for different assumptions. Squares: VER- 
ITAS data. Shaded area: Fermi 68% confidence band. Cir- 
cles: Fermi spectrum. Panel a: Total fluxes for the source 
with unlimited livetime. Panel b: Total fluxes for the source 
with 3-year livetime. 



3 Model Application and EGMF Constraint 

As an example we consider the high-frequency-peaked BL 
Lacertae object (HBL) RGB J0710H-591 located at red- 
shift z — 0.125. The predictions of the total flux as a 
sum of both the direct and cascade photons within the 
instrument point-spread function (PSF) normalized to the 
observed data are shown in Fig. |2] for different EGMF 
strengths and assumptions on source livetime, with a = 
1.5, r = 1/^/1-/32 = 10, and = 25 TeV in Eq.[8] 
The VHE data are from VERITAS measurements [25| 
and the HE data points are extracted from public Fermi 
Large Area Telescope (LAT) data between August 2008 
and January 201 1 using unbinned likelihood analysis in the 
Fermi Science Tools v9rl8p6 with the instrument response 
functions (IRFs) P6_V3_DIFFUS^1 galactic diffuse emis- 
sion model gllJem_v02 and isotropic background model 
isotropicJem_vOfl The ~ 3-year period with simultane- 
ous HE- VHE data sets a lower limit on the livetime of this 
source. 

By requiring that the HE total flux not exceed the Fermi 
LAT measured spectrum, we can roughly see the EGMF 
strength B has a lower limit between 10^^® and 10^^^ 
Gauss for the unlimited livetime case, or between 10^^* 
and 10^^^ Gauss for the livetime assumption of 3 years. A 
more systematic lower limit could be derived by fitting the 
total flux to the measured data points with free parameters 
of normalization /o, index a and cutoff energy i?o. To take 
into account of complexities in blazar modeling (e.g., Il26l ) 
we also fit with a broken power law at 80 GeV and add one 
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Figure 3: Best-fit versus EGMF strength B at different 
livetimes for RGB J0710H-591. 90% and 95% lower limits 
on B are indicated for the 1-year and unlimited livetime 
cases but omitted at other livetimes for clarity. 



more free parameter abreak below the break energy, only to 
find that the resulting lower limit is not greatly affected. 
Restricting the range of a and Ofbreak to be no harder than 
the physically motivated hardness limit 1.5 [27, 281 and the 
range of Eq to be (0. 1 TeV, 100 TeV], we plot the minimum 
from the fit as a function of EGMF strength in Fig. [3] at 
various source livetime limits. As expected all the curves 
converge at low EGMF strengths or large livetimes. 

The EGMF lower limits at different confidence levels are 
derived by finding the point where exceeds its minimum 
value in each curve by Ax^ in Fig. [3] which is just a vari- 
ant of the profile likelihood method for determining confi- 
dence intervals. Two sample confidence levels (90% and 
95%) are given by requiring Ax^ to be 2.72 and 3.84 ||29]| . 
respectively. We show these lower limits versus the blazar 
livetime in Fig. |4] At livetimes below ^ 10^'' years, i.e., 
when the AT constraint is dominating over the Fermi LAT 
PSF constraint on 6c, we have the EGMF lower limit scal- 
ing with AT as i? ^ V AT, consistent with Eqs. [T] and 
|4] The nominal lower limit at 95% confidence level is 
B >2 X 10^^^ Gauss if the source has unlimited livetime 
and B > 3 X 10^^^ Gauss if the source has the minimum 
livetime ~ 3 years. 

4 Conclusions 

The 95% lower limits on EGMF strength inferred by the 
semi-analytic cascade model for both source livetime as- 
sumptions are consistent with the results from full Monte 
Carlo simulations of the cascade in 121] on the same blazar, 

1. Near the completion of this proceeding the Fermi team re- 
leased an updated IRF P6_V1 IJDIFFUSE. We re-analyzed the 
LAT data with the new IRF and applied the cascade model within 
the same analysis framework to find that the final constraint on 
EGMF strength was not significantly affected. 

2. http://fermi.gsfc.nasa.gov/ssc/ 
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Figure 4: EGMF lower limit as a function of blazar live- 
time for two different confidence levels derived from RGB 
J0710H-591 data. Solid black lines: constant B^/Tlines 
agreeing with the curves for small livetimes. 



demonstrating that the model correctly takes into account 
the important geometrical and physical aspects involved in 
the cascade calculation. The limits are conservative as we 
choose a comparably transparent EBL model and also as- 
sume the electron-positron pairs to be always going in one 
coherent magnetic field domain, which is only valid when 
the EGMF coherence length As > 1 Mpc. If As < 1 
Mpc instead, the charged pairs would random walk through 
the EGMF domains |f30l|, and the deflection angle would 
become smaller, causing a larger cascade flux and a more 
stringent lower limit on B. 

Because we ignore cosmological expansion and evolution, 
the model is only applicable to sources within a redshift of 
0.2. However most of the TeV blazars already detected are 
located well in this ranged so we have a pool of candidate 
sources with moderate size. We choose RGB J0710H-591 
as an example because it has simultaneous data in GeV and 
TeV and does not have any variability detected in either en- 
ergy band. The caveat with this one-source study is that the 
EGMF along the line of sight to the source may not be rep- 
resentative of the overall field strength, e.g., when there is a 
filament of intra-cluster magnetic field along the direction. 
Therefore a statistically more reliable constraint should be 
obtained by studying an unbiased set of blazars and using 
the systematic framework presented here to combine the 
results on single sources. As the gamma-ray telescopes 
continue to monitor the sky, more and more blazars with 
simultaneous GeV- TeV baseline spectra will be observed 
and we will be approaching the goal of probing the all-sky 
EGMF 
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